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Abstract. In this paper, the synchronizability problem of dynamical networks 
is addressed, where better synchronizability means that the network synchronizes 
faster with lower-overshoot. The L2 norm of the error vector e is taken as a perfor- 
mance index to measure this kind of synchronizability. For the equilibrium synchro- 
nization case, it is shown that there is a close relationship between the L2 norm of 
the error vector e and the Hi norm of the transfer function G of the linearized net- 
work about the equilibrium point. Consequently, the effect of the network coupling 
topology on the Hi norm of the transfer function G is analyzed. Finally, an optimal 
controller is designed, according to the so-called LQR problem in modern control 
£f) , theory, which can drive the whole network to its equilibrium point and meanwhile 
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minimize the Li norm of the output of the linearized network. 
Keywords, network, synchronizability, L2 norm, H2 norm, LQR problem, 



optimal control. 



1 Introduction 

The topic of synchronizability of dynamical networks has attracted increasing interest recently 
(see [H El El EH EU 111 E] and references therein) , which is mostly referred to as how easy it is for 
the network to synchronize. Technically, this is a problem whether there is a wide range of coupling 
strength in which the synchronization is stable, and the wider the range, the better the synchro- 
nizability of the network. Accordingly, for some cases when the synchronized region is unbounded, 
a better synchronizability means that the synchronization can be achieved with a smaller coupling 
strength [10J. Generally, the synchronizability of a network depends on the underlying coupling 
topology of the network. Studies show that the ability for a network to synchronize is related to 
the spectral properties of the outer coupling matrix of the network [1U|, [TT] and thus influenced 
by the structural properties of the network, such as average distance [H], degree homogeneity [6], 
clustering coefficient [12], degree correlation p], betweenness centrality [3], etc. 
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In this paper, the synchronizability problem of dynamical networks is considered from a different 
point of view: suppose that the coupling strength of a network belongs to the range in which the 
synchronization is stable; then, how fast the synchronization will be achieved? This problem is 
important for the reason that in practical engineering implementation (such as communications via 
chaotic synchronization), synchronization is expected to be achieved not only easily but also swiftly. 
In [7], it is qualitatively pointed out that networks with diagonalizable outer coupling matrices 
may synchronize faster than the ones with non-diagonalizable outer coupling matrices. Clearly, a 
measure is needed for a quantitative description of the swiftness of network synchronization. To 
meet this objective, the L2 norm of the error vector e, denoted as ||e[|2, is taken in this paper as 
a performance index of this kind of network synchronizability. As will be seen later, the quantity 
||e||2 presents a suitable measure of both swiftness and overshoot (referring to the largest difference 
among various node dynamics before the synchronization is achieved): the smaller the quantity 
||e||2, the faster with smaller overshoot the network synchronization. 

Furthermore, as shown by the numerical examples given below, the quantity ||e||2 is influenced 
by the coupling topology of the network. Thus, the investigation on the relationship between ||e||2 
and the network structure is of significance. In this paper, for the case that the synchronous state 
is an equilibrium point, it is pointed out that ||e||2 is upper-bounded by the product of the vector 
2-norm of the initial error vector eo and the H2 norm of the transfer function G(s), denoted as 
||G(s)||2 or simply ||G||2, of the linearized network about the equilibrium point. Thus, the smaller 
the ||G||2, the smaller the ||e||2 as well. As pointed out in [2], the relationship between 1 1 1 1 2 and 
the network structure is quite complicated. Under some assumptions, it is proved in this paper (see 
Theorem 1 and Example 4) that ||G||2 will not increase as the real eigenvalues of the symmetrical 
outer coupling matrix increase. 

For a linear time-invariant system, the linear quadratic regulator problem, or simply the LQR 
problem, is a classical problem in modern control theory. The objective of the LQR problem is 
to find an optimal control law u(t) such that the state x(t) is driven into a (small) neighborhood 
of the origin while minimizing a quadratic performance (L2 performance) index on u and x. In 
fact, the LQR problem is posed traditionally as the minimization problem of the L2 norm of the 
regulator output of the system. In this paper, based on the techniques of the LQR problem, an 
optimal controller design is developed so as to drive the network dynamics onto some homogenous 
stationary states while minimizing the L2 norm of the output of the linearized network. 

The rest of the paper is organized as follows. In Section 2, some preliminary definitions and 
lemmas necessary for successive development are presented. In Section 3, some numerical examples 
are provided to illustrate that the quantity ||e||2 presents a suitable measure of both swiftness and 
overshoot of the network synchronization. For the equilibrium synchronization case, the relationship 
between ||e||2 and the network structure is investigated in Section 4. Based on the results of the 
LQR problem, an LQR optimal controller is proposed in Section 5. The paper is concluded by the 
last Section. 



2 Preliminaries 



Li2[a, b] is an infinite-dimensional Hilbert space, which consists of all square-integrable and 
Lebesgue measurable functions defined on an interval [a, b] with the scalar inner product 

(/,<?>= I f(t)*g{t)dt, 



while if the functions are vector or matrix- valued, the inner product is defined as 

(f,g)= f trace[f(t)*g(t)]dt, 

J a 

where * denotes complex conjugate transpose, and the induced norm is defined as 

ll/lb = VUT), (i) 

for /, g € L 2 [a,b]. 

Consider a continuous-time linear system, 

x = Ax + Bu, y = Cx, (2) 

where x, u and y are the state, input and output of the system, respectively, and A € R nxn , 
B € M. nxm and C G M. lxn are given constant matrices. The transfer function from u to y is 
G(s) = C(sl — A)~ 1 B. If A is stable, then the H2 norm of system ([2]) is represented by the H2 
norm of the transfer function G(s), which is defined by 



l|G( S )|| 2 = \hr I tr^e{G*{ju)G{ju)}dw. 



, 27r 

It can be proved that ||G(s)||2 equals the overall output energy of the system response to the 
impulse input. For computing ||G(s)||, the following formula is convenient. 

Lemma 1 [T3]: If A is stable, then ||G(s)||| = trace(-B T yi?), where matrix Y is the solution to 
the following Lyapunov equation: 

YA + A T Y + C T C = 0. (3) 

Equivalently, 

||G(s)||! = mf {trace(5 T PS) : PA + A T P + C T C < 0} . (4) 

The so-called linear quadratic regulator {LQR) problem is an optimal control problem with a 
quadratic performance {H2 norm) criterion. For the linear time-invariant system 

x = Ax + Bu, x(to) = xq, 

where xq is arbitrarily given, the regulator problem refers to finding a control function u(t) defined 
on [to, T], which can be a function of the state x(t), such that the state x(t) is driven into a (small) 
neighborhood of origin at time T, T < 00. Since every physical system has energy limitation, and 
large control action (even if it is realizable) can easily drive the system out of its valid region, 
certain limitations have to be imposed on the control in practical engineering implementation. For 
these reasons, the regulator problem is usually posed as an optimal control problem with a certain 
combined performance index on u and x. Focusing on the infinite time regulator problem (i.e., 
T — * 00) and without loss of generality assuming to = 0, the LQR problem is formulated as follows: 
Find a control u(t) defined on [0, 00) such that the state x(t) is driven to the origin at t — ► 00 and 



the following performance index is minimized: 
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for some Q = Q*, S and R = R* > 0. Here, R > emphasizes that the control energy has to be 
finite, i.e., u(t) G ^[0, oo). Moreover, it is assumed that 



Then, (ED can be factored as 
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and (J5J) can be written as 
Traditionally, the LQR problem is posed as the following minimization problem: 



min IICsc + -Ditllo- 

u(t)€L 2 [0, oo) 



min 1 1 ?/ 1 1 § 
u(t)ei 2 [o, oo) 



subject to: 



(A4) 



has full column rank for all u. 



x = Ax + Bu, x(to) = xo 
y = Cx + Du. 

For the above LQR problem, the following lemma is useful. 
Lemma 2 [13]: Suppose that in (J7|): 
(Al) (A, B) is stabilizable; 

(A2) D has full column rank with [D D±] being unitary; 
(A3) (C, A) is detectable; 
A - jul B 
C D 

Then, there exists a unique optimal control u = Fx for the LQR problem (f7j). where 

F = -(B*X + DC) 
and X is the stabilizing solution to the following Riccati equation: 

(A - BD*C)*X + X(A - BD*C) - XBB*X + C*D ± D* X C = 0. 
Moreover, the minimized L2 norm of the output y(t) is given by 

|G c x | 

where G c is the transfer function of the system 

f x = (A + BF)x + Ix S(t), z(0_) = 

{ y = (C + DF)x, 
with / being the identity matrix and 5(t) the impulse function. 



(7) 



min 1 1 2/ 1 1 2 

«(t)ei 2 [o, 00) 



1 2 

l'2' 



(8) 
(9) 
(10) 



3 L2 norm performance index of network synchronization 



Consider a network of N identical dynamical nodes, described by 

N 

ii = f(xi) - a^mijFxj, i = 1,2, ••• ,N, (11) 

3=1 

where f(x) governs the dynamics of each individual node, a is the coupling strength, T is the inner 
linking matrix, and M = [rriij] is the outer coupling matrix. 

Let 

e i (t)=x i (t)-x 1 {t), i = !,-■■ ,N, (12) 

and 

e(t) = [ef 4 ••• e T N \ T (13) 
denote the error vector of network (jlip . Then, network (jll|) is said to achieve (asymptotical) 
synchronization if 

e(t) — > 0, as t — > oo. 

Let ||e(t)||2([0,T]), or simply ||e(i)||2 when no confusion may be caused, denote the L2 norm of 
e(t) on a given interval [0, T]. Then, according to (H|), 



(■T 



i/(e,e) = \l J trace[e(i)*e(i)]di. (14) 

In fact, ||e(i)||2 represents the energy of the error signal e on the interval [0, T], which is in 
proportion to the area between the error function e and the time axis. Hence, ||e(t)||2 can be used 
as a quantitative measure of the swiftness and overshoot of the network synchronization. In what 
follows, two examples are first given for illustration. 

Example 1: Suppose that each single node in network (jlip is a Chua's oscillator. In the 
dimensionless form, Chua's oscillator is described by 

x\ = a(-x 1 +x 2 - f{xi)), 

±2 = Xi - X 2 + X3, (15) 
x 3 = -(3x 2 - 7X 3 , 
where /(•) is a piecewise linear function: 

f(xi) = mixi + -(m 2 - mi)(\xi + 1| - \xi - 1|). (16) 

Take parameters a = 9, = 14, 7 = 0.01, mi = —0.714, and m 2 = —1.14, so that Chua's 
oscillator (fl~5|) generates a double-scroll chaotic attractor. Set T = diag(l, 1, 1) in (fTTI) . Fig. 1 
shows the different synchronization performances of network (jlip with the same coupling strength 
a [a = 6) but different coupling configurations. The corresponding values of ||e(t)||2 are computed 
numerically as given in Table 1. 



Nearest-neighbor coupling 



Star-shaped coupling 



Global coupling 





Nearest-neighbor coupling Star-shaped coupling Global coupling 

Fig. 1 Synchronization of the state variables of network (TTTj) with 20 nodes and 
the same coupling strength in different coupling configurations. 
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Nearest-neighbor 


Star-shaped 


Global 


1Mb 


9.2744 


0.8030 


0.4142 



Table 1 Values of II ell 



As to the case of equilibrium synchronization of network (jlip . since the synchronous state s is 
known, the error vector can be defined in the following way: 

ei = Xi(t) - s, i = 1, ■ • • ,N, 
e(t) = [e\e\ ■■■ e T N ] T . 

Example 2: Consider a Lur'e system, 

f x x = (A x - 2B 1 C 1 )x 1 + Bi/ifc/i), 
\ yi = Cixi, 
where X\ is the state, y\ is the measured output, 
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and the nonlinear function /i(yi) = \yx + 1| — |yi — 1|. 

A network with system (|18p as individual nodes is given as follows: 
f x = (I N ® (Ai - 2B!Ci) - oM ® T)x + (Jjv 
1 V = {1n® Ci)x, 



(17) 



(18) 



(19) 



where x = (x{ , • • • ,x^) T , y = (y x , ■ ■ ■ ,y N ) T and f(y) = (/i(yi), • • • , /tv(jmO) T - 

Network (fT9j) can be viewed as a large-scale system with measured output and feedback. Assume 
that system (|19p is observable. Then, the states of network nodes achieve synchronization if and 
only if all the outputs of network nodes achieve synchronization. Thus, one only needs to consider 
the outputs of network (fT9j) . 

Replace X{ with yi and set s = in (fT7|) . Then, a = yi and e = y. Let a = 10, 6 = 3, 61 = 
0, 62 = 1, c\ = C2 = 1, and cr = 2 in (|18|) . Fig. 2 shows the different performances of outputs of 
network (|19p with the same coupling strength a (cr=2) but different coupling configurations. The 
corresponding values of ||y||2 are listed in Table 2. 






Nearest-neighbor coupling Star-shaped coupling Global coupling 

Fig. 2 Synchronization of the outputs of network (fT!))) with 20 nodes having 
the same coupling strength but different coupling configurations. 



Coupling 


Nearest-neighbor 


Star-shaped 


Global 


Wvh 


0.7296 


0.5768 


0.3990 



Table 2 Values of ||y|| 2 . 



□ 

Remark 1: As indicated by the above examples, in the same situation of stable synchronization, 
the performances of synchronization of networks can be quite different. 

Remark 2: It is also clear that swiftness and overshoot are important indexes for describing the 
synchronous behaviors. The L2 norm of the error vector e, i.e., ||e||2 as defined in (|13|) or (|17p for the 
equilibrium synchronization case, can properly measure the swiftness and overshoot of the network 
synchronization: the smaller the ||e||2, the faster with lower overshoot the network synchronization. 
Thus, the quantity ||e||2 can be taken as a performance index of network synchronizability. 

Remark 3: The quantity ||e||2 is influenced by many factors of the network, such as 

(1) network structure, particularly the coupling strength a and the outer coupling matrix M; 

(2) dynamical components, particularly the individual dynamics determined by /, the syn- 
chronous state s(t), and the inner linking matrix T. 



4 Local synchronization to network equilibrium point 



In the literature where the eigenvalues of the outer coupling matrix are used to measure the 
network synchronizability, it is a topic of great interest to investigate the relationship between the 
eigenvalues and the network structure thereby finding suitable ways to enhance the synchroniz- 
ability. In this section, for the case of equilibrium synchronization, the relationship between the 
quantity ||e||2 and some network parameters is explored. 



4.1 HGH2 as a constraint of ||e|| 2 

Suppose that each single node in network (fTT]) has a measured output y = Cx. Then, the 
equations of network (jlip can be written as follows: 



N 

Xi = f(xi) - a Yl mijTxj, 2^(0) 
3=1 

Hi = Cxi, 



XiO, 



1. 



,N. 



(20) 



If C = I n , then yi = Xi is just the state of the ith node. 

Let s denote an equilibrium point of the individual node, satisfying f(s) = 0. Then, the linearized 

equations of (|2U|) about the synchronous solutions Xi — S 3X6 3S follows: 

N 

Vi = D f( s )Vi ~ a E mijTrij, ^(0) = x i0 - s = rn , . 

3=1 i = l,"' ,N, (21) 

&i = Crji, 

where rji (rji = Xi — s) and are the state error vector and the output error vector to the ith node, 
respectively. Viewing the impulse function r]io6(i) as an input to system (f2"T|) . system (f2"T|) can be 
equivalently written as 

( N 

i]i = Df(s)rji -ctYj raijTrij + i] i0 5(t), 7ft (0_) = 0, 

j=i 1 -- 

ei = Crji, 

Using the Kronecker product, the error system (|22p can be rewritten as 

f r) = (I N ® Df(s) -oM® F)rj + {I N (g> J„)«(t), 
1 e = (Jjv (8 C)r?, 



(22) 



(23) 



where 7/ = [r]J 77J • • • ??^] T , e = [ef e T 



Q TlT 



and «(i) = r]oS(t). 



As in Example 2, suppose that system (|23|) is observable. Then, ||e||2 can be taken as a measure 
of the swiftness and overshoot of the synchronization of network (|20p . Let G(s) denote the transfer 
function of the error equation (|23p from u to e. Then 

G(s) = (I N ® C)[s(/jv ® / n ) - (Jjv ® - <rM ® T)]- 1 ^ ® In). (24) 

Let the transfer function G(s) be given in (|24p . then Lemma 3: The inequality 

||e|| 2 < IIGH2IMI2 (25) 
holds, where ||e||2 is the L2 norm of the error vector e(t) on the interval [0, 00), ||Gr||2 is the H2 
norm of the transfer function G(s), and 1 1 970 1 1 2 is the vector 2-norm of the initial error vector t]q. 

Proof: Since u(t) = r]o5(t), e(t) = g(t)rjo, where g(t) denotes the corresponding bilateral Laplace 



transform of G(s). Then, by Parseval's identity, 

lleWlb = \\g{t)voh = \\Gmh ^ 1 1 <^ 1 1 2 1 1 »7o 1 1 2 - 



□ 



Thus, the quantity ||e||2 is upper-bounded by the product HGI^H^olh- Since 770 is the given initial 
error vector, ||G||2 can be taken as a constraint of ||e||2- In fact, as introduced in Sec.l, ||G||2 equals 
the overall output energy of the system response to the impulse input. 

Remark 4: The advantages of using the quantity ||G||2 include: 

(1) ||G||2 can be numerically computed; 

(2) the synchronizability is affected by many factors of a network, while \\G\\2 can be seen as an 
overall reflection of these network factors consisting of both structural and dynamical ones. 

Example 3: Consider Example 2 again. The linearized equation of network (|19[) about the 
equilibrium point s = [0 0] is given as follows: 

f 7) = (Jjv ® A x - oM ® V)r\ + (I N ® I n )voS(t), v(P-) = 0, 
1 e= {I N ®C)r]. 

The corresponding values of ||G||2 with the three different network configurations are listed in 
Table 3. 



Coupling 


Nearest-neighbor 


Star-shaped 


Global 


l|G||2 


2.5091 


2.4383 


1.2565 


ll e l|2 


0.7296 


0.5768 


0.3990 



Tabic 3 Values of ||G||2 of network (Tl9|) with different network configurations. 

4.2 Relationship between ||G||2 and network structure 

Example 3 shows that ||G||2 is influenced by the network configurations. Thus, the relationship 
between ||G||2 and the network structure is a problem deserving further investigation. 

In this section, it is always assumed that the outer coupling matrix M is symmetrical. Then, 
there exists a unitary matrix U G R JVxAr such that M = £/At/ _1 , where A = diag(Ai, A2, • • • , \n) 
is a diagonal matrix with the diagonal entries A«, i = 1, • • • , N, being the eigenvalues of matrix 
M. Let r\ = (U <S> I n )£, u = (U ® I n )u>, and z = (U^ 1 <£> I{)e. Then, the error equation (|23|) can be 
rewritten as follows: 

{ i = {I N ® Df(s) - aA ® r)f + (I N ® I n )u, 
\ z=(I N ® G)^. 
Note that system (p7|) is composed of N uncoupled subsystems: 

f it = (Df(s) - aXiT% + I n cj h 
I Zi = C£i, 



,N, 



(27) 



(28) 



where loi is the ith component of the input vector u> = (U 1 <£> I n )u = (U 1 <8> I n )rjo5(t). Let 

Ek = diag(0 5 • • • ,0,1,0,-. ■ ,0), 

be a matrix with the ith diagonal entry being 1 and all the other entries being zero. Then, 
oji = (Ei (g> I n ){U- 1 I n )uj = (EiU- 1 ® I n )VoS(t) and uj i0 = {Eft- 1 ® I n )r] . 



From (I28p . the condition for ensuring the stable synchronous state 

Xi = X2 = • ■ ■ = Xn = S 

is that the ./V matrices 

Df(s)-aXiT, i = l,.-.,N, (29) 

are all stable. 

Let T(s) denote the transfer function of (|27p from uj to z. Then 

T(s) = (I N ® C)[s(I N ® / n ) - (Jjv ® D/(s) - aA ® T)]- 1 ^ (8 J n ) 
= (C/- 1 ^/ / )G(s)([/^/ m ). 
Since both (t/ -1 ® //) and (£/ ® I m ) are unitary transformations, one has 

llrOOIIa = ||G( S )|| 2 . (30) 

Let Tj(s) = C[sJ n — (Df(s) — <rAjr)] _1 / n denote the transfer function of the ith. subsystem in ([28]). 
Then 



and 

iv 



r(«)iii = Eii r *wiii- ( 31 ) 

Assumption 1: Suppose that the outer coupling matrix M is symmetrical, diffusive and irre- 
ducible, with the off-diagonal entries < and the diagonal entries ma = ~Ylf=i j^i m ij) f° r 
i, i = 1,2, • • • ,N. 

The outer coupling matrix M satisfying Assumption 1 has a zero eigenvalue of multiplicity 1, 
and its other eigenvalues are all positive. 

Theorem 1: Suppose that the inner linking matrix T = kl, where k > is a constant and I 
the identity matrix. Suppose that the outer coupling matrix M satisfies Assumption 1 with the 
eigenvalues given as follows: 

= Ai < A 2 < ■ •• < Atv. (32) 

Then 

\\Ti(s)\\l< 11^)111, (33) 

where 1 < j < i < N, and 

ITOIIl > \\Tl(s)\\l + (N- l)\\T N (s)\\l > N\\T N (s)\\l 
\\T(s)f 2 < HTiWHl + {N- 1)||T 2 ( S )||| < ^||r x ( fl )|||. 
Proof: Let P be an arbitrary positive definite matrix such that 

P(Df(s) - aXjT) + (Df(s) - aXjTfP + C T C < 0. 
By the assumption that T = kl, one has 

PT + T T P > 0. 



(34) 



□ 



Then 

P(Df(s) - aX l T) + (Df(s) - aXiT) T P + C T C 

= P(Df(s) - aXjT) + (Df(s) - aXjTfP + C T C - a(Xi - Xj)(PT + T T P) 
< -a{Xi - Xj){PT + T T P) < 0. 
The above inequality and Lemma 1 together leads to the assertions of the theorem. 

Remark 5: Since Ai = 0, ||Ti(s)||2 in (|34f) represents the H2 norm of the transfer function of 
each individual node in the network. 

Remark 6: Theorem 1 provides a simple way for comparing ||G||2 of a network with its different 
possible structures. It can be deduced from Theorem 1 that ||G||2 will not increase as the eigenvalues 
of the outer coupling matrix increase. 

Example 4: Consider Example 2 again. The eigenvalues of the global coupling matrix M g i Q are 
Ai = and A2 = ■ • • = Xn = N, while the eigenvalues of the star-shaped coupling matrix M sta are 
Ai = 0, A2 = • • • = Aat_i = 1 and A at = N. Let G g i Q and G s t a denote the transfer functions of the 
global coupling and the star-shaped coupling networks, respectively. Then, by Theorem 1, 

= \\T l {s)\\l + {N-l)\\T x=N \\l 

<l|r 1 (-)||l + ||2wil3 + (^-2)||2WHl 



\C II 2 
\ K - T glo\\2 



\T l|2 
1 1 glo\\2 



\T II 2 
I x sta || 2 



\r II 2 

\^ sta ||2- 



It is consistent with the numerical results given in Example 3. 

For Example 2, through integral computations, an analytical expression of ||Tj||2 can be obtained 
as follows: 



2(cf + 4)a 2 Xf + [(3cf + 4)b + 2cic 2 (l - a)]<rA; + (ab - c 2 a) 2 + (cf + 4)a + c( 



(35) 



" 112 4a 3 Af + 6ba 2 Xf + (4a + 2b 2 )aXi + 2ab 

By (j35|) and (f3T|) . Table 4 gives the different values of ||G 9 ; ||2 and ||G S 4 a |j2 of network (fT9]) as the 
node number N increases. 
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10 


100 


1000 


00 


\\Gglo 2 


1.0801 


1.3190 


1.4492 


1.4696 


1.2910 


\\G s ta 2 


1.0801 


2.2619 


6.9840 


22.0434 


00 



Tabic 4 Values of ||G||2 of network (fT9|) with different network configurations. 



□ 



5 Optimal controller design 

So far, the pinning control strategy is extensively used for achieving synchronization of dynamical 
networks [U [8]. The main advantage of the pinning control strategy is that only a few network nodes 
are needed to be controlled. AS far as the control effects (referring to the swiftness and overshoot 
of the synchronization) and the control cost (represented by the L2 norm of the control input) are 
concerned, however, pinning control may not be the best choice. It is revealed by the LQR problem, 
as illuminated in ([7]), when matrices C and D are properly selected, both the control effects and 
the control cost can be simultaneously evaluated by the L2 norm of the measured output y of the 



controlled network. The smaller the ||y||, the better the control effects and the lower the control 
cost. In this section, based on the optimal solution to the LQR problem, an LQR optimal controller 
is designed for network (jlip . which can drive the network onto some homogenous stationary states 
while minimizing the quantity ||y||. 

Suppose that the controlled network is given as follows: 

N 

Xi = f(xi) - a Yl rriijTxj + Bu { 

0=1 i = !,■■■ ,N, (36) 

Ui = Cxi + Duj, 

where B e M. nxm , C G R lxn , D E M. lxm are given constant matrices, u = (uj , • • • ,uJj) T is the 
controller to be designed, and y = (yf , • • • , yJj) T is the measured output of the LQR problem. As 
in Sec. 4, let s denote an equilibrium point of the individual node and let r\ = (r]f , • • • , r]J r ) T , where 
r]i = Xi — s, i = 1, • • ■ , N, be the state error vector. 

Suppose that the controller is a feedback law of the state error, i.e., u = Fij, where F is 
the feedback gain matrix to be determined. Then, by the Kroneck product, the equation of the 
linearized system of (|36p about the synchronous solution x\ = ■ ■ ■ = xn = s is given as follows: 

{77 = A c r\ + B c u, , s 

y = C c r) + D c u, 

where A c = I N ® Df(s) - crM <g> T, B c = I N ® B, C c = I N C and D c = I N ® D. Furthermore, 
suppose that the matrices A c , B c , C c and D c satisfy the assumptions (A1-A4) as given in Lemma 
2. Then, by Lemma 2, the feedback gain F for the LQR problem of the linearized system (|37p is 
obtained as 

F = —(B*X + D C C C ), (38) 
where X is the stabilizing solution to the following Riccati equation: 

(A c - B C D*C C )*X + X(A C - B C D* C C C ) - XB C B*X + C*{D C ) ± {D C )* X C C = 0. (39) 



Remark 7: The optimal controller u = Frj with F given in (|38p has the property that it controls 
network (|36p to the synchronous state x\ = ■ ■ ■ = xn = s and meanwhile minimizes the L2 norm 
of the output y of the linearized system (|3"Tj) . 

As the node number iV of network (|36|) increases, the direct computation of the feedback gain 
F as given in (|38p will become harder. For the case when the outer matrix M is symmetrical, the 
feedback gain F can be alternatively given in terms of the iV uncoupled subsystems. The main 
advantage of this approach is that the decentralized feedback gain Fi will be given only based on 
the information of the ith subsystem, i = 1, 2, • • ■ ,N. 

By using the unitary transformation, as used in Sec. 4, system (|37p becomes 

f £ = (I N ® Df{s) - a A ® T)£ + (Jat ® B)a;, 
\ J8 = (JiV ® C)^ + (JiV ® 
where r\ = (U®I n )£, u = (U ® I n )u, and 2; = (f/ -1 ®I/)y. The uncoupled controlled subsystems 
are as follows: 

^ ^ i = l,..-,JV, (41) 



where A4 = Df(s) 



a\iT. By the analysis in Sec. 4, one has 

N 

II ||2 _ || || 2 _ || i|2 

\\y\\2 ~ \\ z \\2 — 11^112- 

i=l 



Suppose that the constant matrices Ai, B, C and D in (|41[) satisfy the assumptions (A1-A4) as 
given in Lemma 2, for i = 1, • • • , N. Then, the controller to the ith subsystem can be given as 
follows: 

Fi :=-(B*X l + DC), (42) 
where Xi is the stabilizing solution to the following Riccati equation: 

(Df(s) - XiT - BD*C)*X l + Xi(A- XiT - BD*C) - XiBB*X t + C*D ± D* ± C = 0, (43) 
for i = 1, 2, • • • ,N. Moreover, the feedback gain matrix F is given by 

Fx 



F N 



(U* I n 



(44) 



Example 5: Consider Example 1 again. Suppose that the network is with the global coupling 
configuration. The objective is to design a controller such that network (jlip with Chua's oscillators 
can be driven to the synchronous solution x\ = • • • = xn = s = [0 0] T . 

It is easy to deduce that 

—a(m2 + 1) a 
Df(s)= 1 -1 1 

-13 - 7 

Consequently, the linearized equation of the controlled network is given by 

J fj=(I N ® Df(s) -aM® T)rj + (I N ® B)u, 
\ V = (In® C)r] + (Jjv (8 D)u, 
where the constant matrices B = C = D = 1% in this example. 

By the above analysis, the LQR optimal controller can be designed by solving the Riccati equation 
(|39p . For comparison, a pinning control strategy is also considered, where 6 nodes are randomly 
selected to be pinned by controllers it, = —a ■ d ■ I^rji with a = 1 and d = 20. 

Fig. 3 shows the different performances of the synchronous behavior of the two controlled net- 
works, by the LQR optimal controller and the pinning controller, respectively. The corresponding 
values of \\y\\2 are listed in Table 5. 



(45) 





LQR optimal control pinning control 

Fig. 3 Different effects of the LQR and pinning controllers. 



Controller 


LQR optimal control 


Pinning control 


Wvh 


0.1309 


4.9799 



Table 5 Values of ||y|| 2 . 



□ 



6 Conclusion 

In this paper, synchronizability of dynamical networks is considered based on some new measures: 
the swiftness and overshoot of the network synchronization. The quantity ||e||2, which represents 
the Li norm of the synchronization error vector e(t), is taken as the performance index of this kind 
of synchronizability. It has been shown by several numerical examples, ||e||2 presents a suitable 
measure of both swiftness and overshoot of network synchronization: the smaller the values of ||e||2, 
the faster with smaller overshoot the network synchronization. For the case when the synchronous 
state is an equilibrium point, ||e||2 is upper-bounded by the product of the vector 2-norm of the 
initial error vector eo and the Hi norm of the transfer function G(s), denoted as ||G(s)||2, of the 
linearized network about the equilibrium point. The relationship between ||G(s)||2 and the network 
structure has also been discussed. Under some assumptions, it has been proved that 1 1 1 1 2 will 
not increase as the real eigenvalues of the outer coupling matrix increase. Finally, based on the 
techniques of the LQR control theory, an optimal controller has been suggested to drive the network 
onto some homogenous stationary states, which , in the mean time, can minimize the Li norm of the 
output of the linearized network. Further research along this direction seems to be quite promising 
as long as the network energy and performance are concerned, therefore deserves further efforts. 
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